Three-body Recombination of ^Li Atoms with 
Large Negative Scattering Lengths 
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The 3-body recombination rate at threshold for distinguishable atoms with large negative pair 
scattering lengths is calculated in the zero-range approximation. The only parameters in this limit 
are the 3 scattering lengths and the Efimov parameter, which can be complex valued. We provide 
semi-analytic expressions for the cases of 2 or 3 equal scattering lengths and we obtain numerical 
results for the general case of 3 different scattering lengths. Our general result is applied to the three 
lowest hyperfine states of ®Li atoms. Comparisons with recent experiments provide indications of 
loss features associated with Efimov trimers near the 3-atom threshold. 
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Atomic gases allow the experimental study of super- 
fluidity in systems in which the fundamental interactions 
are simple and experimentally controllable. In the case 
of fermionic atoms with two spin states, there have been 
extensive investigations of the crossover from the BCS 
mechanism (Cooper pairing of atoms) to the BEC mecha- 
nism (Bose-Einstein condensation of diatomic molecules) 
Fermionic atoms with three spin states open up the 
possibility of new superfluid phases and new mechanisms 
for superfluidity 0, y, 0, 0] • The first experimental stud- 
ies of such a system were recently carried out using the 
three lowest hyperfine states of ^Li atoms d, 0] • 

Fermionic atoms with three spin states also open up 
new possibilities in few-body physics. For the discussion 
of 3-body observables, they can simply be considered as 
three distinguishable atoms and their fermionic nature 
plays no special role. If the pair scattering lengths are 
large compared to the range of the interactions between 
the atoms, a remarkable set of 3-body phenomena are 
predicted. If the 3 scattering lengths are infinitely large, 
there is an infinite sequence of 3-atom bound states called 
Efimov trimers with a geometric spectrum and an accu- 
mulation point at the 3-atom threshold [8]. Low-energy 
3-body phenomena governed by discrete scale invariance 
are generally referred to as Efimov physics 0, [l3| ■ The 
first experimental evidence for Efimov physics were ex- 
periments with ultracold ^^'^Cs atoms by Grimm and 
coworkers in which they observed dramatic dependence 
of the 3-body recombination rate and the atom-dimer 
relaxation rate on the scattering length 
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In this Letter, we present calculations of the 3-body re- 
combination rate at threshold for distinguishable atoms 
with large negative pair scattering lengths in the zero- 
range limit. We provide semi- analytic expressions for 
the cases of 2 or 3 equal scattering lengths and we obtain 
numerical results for the general case of 3 different scat- 
tering lengths. We apply our general results to the three 



lowest hyperfine states of ®Li atoms and compare with 
recent 3-body recombination rate measurements 0, 0] • 

We consider an atom of mass m with three distinguish- 
able states that we label 1, 2, and 3 and refer to as spin 
states. We denote the scattering length of the pair i and 
j by either = aji or Ok, where (ijk) is a permutation 
of (123). The rate equations for the number densities Ui 
of atoms in the three spin states are 



dt 



-K^nin2n3. 



(1) 



By the optical theorem, the event rate constant ^3 in 
the low-temperature limit can be expressed as twice the 
imaginary part of the forward T-matrix element for 3- 
atom elastic scattering in the limit where the momenta of 
the atoms all go to 0. Using diagrammatic methods, the 
T-matrix element for elastic scattering can be expressed 
as the sum of 9 amplitudes corresponding to the 3 pos- 
sible pairs that are the first to scatter and the 3 possible 
pairs that are the last to scatter. For small collision ener- 
gies, the leading contributions to those amplitudes come 
from the S-wave terms, which we denote by Aij{p,p'), 
where p (p') is the relative momentum between the pair 
that scatters first (last) and the third atom labelled i (j). 
The rate constant in Eq. ([1]) is 



= ^ aiajlinAtj (0, 0), 

i,3 



(2) 



where the sums are over i,j — 1,2,3. The amplitudes 
Aij{p,p') can be calculated in the zero-range limit by 
solving 9 coupled integral equations that are generaliza- 
tions of the Skorniakov-Ter-Martirosian (STM) equation 
[ist . To determine ImAij (0,0), it is sufficient to solve the 



2 



9 coupled STM equations for Aij {p, 0) : 

X / dqQ{q/p)D,{q)Akiq,0), (3) 

where 

Dkiq) = i-l/ak + VSq/2)-\ (5) 

and A is an ultraviolet cutoff. The solutions to Eqs. ([3]) 
are singular as p — > 0. The singular terms, which 
are proportional to 1/p, and Inp, appear only in 

Re Aij {p, 0) for real p and can be derived by iterating 
the integral equations ■ Since Im Aij {p, 0) must be ex- 
trapolated to p = 0, it is useful to transform Eqs. ([3]) into 
coupled STM equations for amplitudes Aij {p, 0) obtained 
by subtracting the singular terms from Aij{p,0). For 
p <C A, the solutions depend log-periodically on A with a 
discrete scaling factor e^/*° w 22.7, where Sq = 1.00624. 
The dependence on the arbitrary cutoff A can be elim- 
inated in favor of a physical 3-body parameter, such as 
the Efimov parameter defined by the spectrum of Efi- 
mov states in the limit where all 3 scattering lengths are 
infinitely large [§]: 

^2./so\ (ai2 = - aai = ±oo). 

/ m 

(6) 

If we restrict A to a range that corresponds to a multi- 
plicative factor of 22.7, then A differs from only by 
a multiplicative numerical constant. Thus we can also 
simply take A as the 3-body parameter. 

If ttij > 0, there is a contribution to K3 from 3-body 
recombination into the shallow dimer whose constituents 
have spins i and j and whose binding energy is 'h?' / [mafj) . 
If ai2, 023, and 031 are all negative, there are no shal- 
low dimers. The solutions Aij{p,Q) to the coupled STM 
equations in Eq. ([3]) are all real-valued in this case, so 
the rate constant if 3 in Eq. ^ is predicted to be 0. 

If there are deeply-bound diatomic molecules {deep 
dimers) in any of the three 2-body channels, there are 
also contributions to K3 from 3-body recombination into 
the deep dimers. If all 3 scattering lengths are negative, 
these are the only contributions to K^. The coupled STM 
equations in Eq. ^ do not take into account contribu- 
tions from deep dimers. The inclusive effect of all the 
deep dimers can be taken into account by analytically 
continuing the Efimov parameter to a complex value 
(li ]: — > K, exp(i77,/so), where 77, is a positive real 
parameter. Making this substitution in Eq. ([S]), we find 
that the Efimov states acquire nonzero decay widths de- 
termined by rj^ . If we use the ultraviolet cutoff A as the 
3-body parameter, the inclusive effects of deep dimers 




FIG. 1: (Color online) The coefficients C scaled by (a23/a)~^ 
(upper curve) and D (lower curve) in Eq. (O as functions of 
a-izja for the case of two equal negative scattering lengths a 
and a third negative scattering length 023. 

can be taken into account by changing the upper limit of 
the integral in Eq. ([3]) to A exp(i7;,/so), so the path of in- 
tegration extends into the complex plane. Having made 
this change, the solutions Aij (p, 0) are complex- valued 
even if ai2, 023 and 031 are all negative. The rate con- 
stant K3 in Eq. (2) is a function of the scattering lengths 
ai2, 023, and 031 and the 3-body parameters A and 77* 
and it vanishes as ry* — > 0. It gives the inclusive rate for 
3-body recombination into all deep dimers. 

We focus our attention on cases in which all scattering 
lengths are negative, so the only recombination channels 
are into deep dimers. We first consider the case of 3 equal 
scattering lengths: 012 — 023 = 013 = a < 0. In this case 
Eq. ([3|) reduces - after summing over i a nd 7 - to the STM 
equation for identical bosons. In Ref. [l5|, Braaten and 
Hammer deduced an analytic expression for the 3-body 
recombination rate constant for identical bosons with a 
large negative scattering length a: 

^ _ 167r^Csinh(27?,) ha^ 

sin^[soln(D|a|K*)] -|-sinh^77* m 

where sq — 1.00624, C, and D are numerical constants. 
This formula exhibits resonant enhancement for a near 
the values (e'^/*'')"(DK*)~^ for which there is an Efimov 
state at the 3-body threshold. Fitting our numerical 
results for K^/a* as functions of aA and 77*, we deter- 
mine the numerical constants to be C = 29.62(1) and 
D = 0.6642(2). These values are more accurate than 
previous results for identical bosons A separate cal- 
culation of the spectrum of Efimov states in the limit 
a — *■ ±00 with 77, = is necessary to determine the re- 
lation between the Efimov parameter and the ultraviolet 
cutoff: = 0.17609(5)A. 

We next consider the case of 2 equal negative scattering 
lengths and a third that vanishes: ai2 — 013 = a < 0, 
0-23 = 0. In this case with only two resonant scattering 
channels, sq = 0.413698 and the discrete scaling factor 
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FIG. 2: (Color online) The scattering lengths for the three 
lowest hyperfine states of ^Li as functions of the magnetic 
field B [10]. The vertical scale changes by a factor of 10 at 
B — 600 G. The two vertical dotted lines mark the boundaries 
of a region in which |ai2| > 2 Ivdw- The three vertical dashed 
lines mark the positions of the Feshbach resonances. 



is e'^^^° ~ 1986. Eq. ([7]) again gives an excellent fit to 
our numerical results and we determine the numerical 
constants as C = 0.8410(6) and D = 0.3169(1). 

We now consider the case of 2 equal negative scattering 
lengths and a third that is unequal: ai2 = ai3 = a < 0, 
023 < 0. Eq. ([7]) with sq — 1.00624 continues to provide 
an excellent fit to our numerical results. The fitted val- 
ues of C and D are shown as functions of 023 /a in Fig. [T] 
For 1 023 1 ^ the coefficients seem to have the limit- 
ing behaviors C w 10.88(2) (023/0)^ and D « 1.30(1). 
Their limiting behaviors for |a23| <C \a\ do not seem to 
be simple power laws. This is not surprising, because the 
discrete scaling factor 22.7 changes to 1986 when 023 — 0. 

Finally we consider the general problem of 3 different 
negative scattering lengths, for which we can obtain nu- 
merical results for given values of 012, 023, and 013. We 
apply our method to ^Li atoms in the three lowest hy- 
perfine states \ f,mf): |1) = ll,-!-^), |2) = |^,-^), and 
|3) = ||,— |). The 3 pair scattering lengths 012, 023, 
and ai3 are shown as functions of the magnetic field in 
Fig. [5] [161]. They have Feshbach resonances near 834 G, 
811 G, and 690 G, respectively [I?!- 



The zero-range 
approximation should be accurate if lai2l, 1023], and 
]ai3] are all much larger than the van der Waals length 
^vdw = {mCQ/ti^y^^ , which is approximately 62.5 oq 
for ^Li. There are two regions of the magnetic field in 
which all 3 scattering lengths are negative and satisfy 
]ay ] > 2^vdw: a low-field region 122 G < i? < 485 G and 
a high-field region B > 834 G. In the low-field region, the 
smallest scattering length is ai2 and it achieves its largest 
value —290 oq = —4.6 ^vdw near 320 G. The zero-range 
approximation may be reasonable near this value of B. 
In the high-field region, the smallest scattering length 
is ai3. It increases from —3285 uq at B — 834 G to 
—2328 ao ~ —37 ^vdw at 1200 G. Thus the zero-range 
approximation should be very accurate in this region. We 
emphasize that the 3-body parameters and 77* need 
not be the same in the two universal regions, since there 
are zeroes of the scattering lengths between them. 
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FIG. 3: (Color online) The 3-body recombination rate con- 
stant K3 as a function of the magnetic field B. The two ver- 
tical dotted lines mark the boundaries of the region in which 
]ai2] > 2 ^vdw. The solid squares and dots are data points 
from Refs. ^ and [J, respectively. The curve is a 2-parameter 
fit to the shape of the data from Ref. @. 



The 3-body recombination rate K3 for ^Li atoms in the 
three lowest hyperfine states has recentlybeen measured 
by Jochim et al. and by O'Hara et al. [7[. Their results 
are shown in Figs. [3] and ID In Ref. if 3 was measured 
for each of the three spin states separately. Those results 
have been averaged to get a single value of K3 at each 
value of B. Both groups observed dramatic variations in 
K3 with B, including a narrow loss feature near 130 G 
and a broader loss feature near 500 G. 

The narrow loss feature and the broad loss feature ob- 
served in Refs. 0] appear near the boundaries of the 
low-field region in which all 3 scattering lengths satisfy 
\aij\ > 2^vdw- The zero-range approximation is question- 
able near the boundaries of this region. We nevertheless 
fit the data for K3 in this region by calculating the 3- 
body recombination rate using the J5-dependence of 012, 
a23, and 013 shown in Fig. [21 while treating A and ry* 
as fitting parameters. Since the systematic error in the 
normalization of K3 was estimated to be 90% in Ref. @ 
and 70% in Ref. [T], we only fit the shape of the data 
and not its normalization. A 2-parameter fit to the data 
from Ref. |B] in the region 122 G < S < 485 G gives 
A = 436 Oq ^ and ry* = 0.11. The fit to the shape of the 
narrow loss feature is excellent as shown in Fig. [3| Hav- 
ing fit the position and width of the loss feature feature, 
the normalization of K3 is determined. In the region of 
the narrow loss feature, the prediction for K3 lies be- 
low the data of Ref. Q by about a factor of 2, which is 
well within the systematic error of 90%. The excellent fit 
to the shape of the narrow loss feature and the predic- 
tion of the normalization of consistent with the data 
suggests that this loss feature may arise from an Efimov 
state near the threshold for atoms in spin states 1,2, and 
3. As shown in Fig. [31 our fit predicts that should 
be almost constant in the middle of the low-field region 
and that there should be another narrow loss feature at 
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FIG. 4: (Color online) The 3-body recombination rate con- 
stant Ki as a function of the magnetic field B. The three 
vertical dashed lines mark the positions of the Feshbach res- 
onances. The solid dots are data points from Ref. 0]. The 
curve is a 2-parameter fit to the last two data points. 

its upper end near 500 G. The data from both groups in 
Fig. [3] increases monotonically in the middle of the low- 
field region and, instead of a narrow loss feature, there 
is a broad loss feature near the upper end of this region. 
We are unable to get a good fit to the slope of log Kj, in 
the middle of the low-field region or to the shape of the 
broad loss feature by adjusting A and 77*. 

In Ref. 0], the 3-body recombination rate was also 
measured at higher values of the magnetic field. They in- 
clude three data points in the region B > 834 G, where 
all 3 scattering lengths are extremely large and nega- 
tive. If the central values of the last two data points 
are used to determine the 3-body parameters, we ob- 
tain A = 37.0 Oq ^ and = 2.9 x 10^"*. As shown 
in Fig. m this fit predicts the resonant enhancement of 
the 3-body recombination rate near 1160 G. If we al- 
low for the systematic error by increasing or decreasing 
both data points by 70%, the position of the resonance 
does not change, but 77* increases to 5 x lO"'^ or de- 
creases to 9 X 10^^, respectively. If we take into account 
the statistical errors by increasing or decreasing the data 
points by one standard deviation, the position of the res- 
onance can be shifted downward to 1109 G or upward 
to 1252 G. Thus it might be worthwhile to search for an 
Efimov resonance in this region. If such a feature were 
observed, measurements of its position and width would 
determine accurately the two 3-body parameters and 
77* . Our equations could then be used to predict the total 
3-body recombination rate in the entire universal region 
B > 610 G, including the regions where 1, 2, or 3 of the 
scattering lengths are positive. Note that the third-to- 
last data point in Fig. [4] shows no sign of the large in- 



crease in K3 near the Feshbach resonance at 834 G that 
is predicted by our fit. However the measurement of K3 
involves a model for the heating of the system, and the 
failure of our fit at 835 G might be attributable to the 
breakdown of that model near the Feshbach resonance. 

In summary, we have calculated the recombination rate 
of three distinguishable atoms with large negative pair 
scattering lengths in the zero-range limit. We have pro- 
vided simple semi-analytical expressions for the rate if 
2 or 3 scattering lengths are equal. Using our general 
result for 3 unequal scattering lengths, we showed that 
the narrow 3-body loss feature for ^Li atoms with three 
spin states [1, 01 may be attributed to an Efimov state 
near threshold. Our results provide a starting point for a 
quantitative understanding of the 3-body loss rates and 
the unambiguous identification of Efimov physics in these 
systems. 
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